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Gamma-ray bursts (GRBs) are one of the candidates of ultra-high-energy (> 10 18 ' 5 eV) cosmic- 
ray (UHECR) sources. We investigate high-energy cosmic-ray acceleration including heavy nuclei in 
GRBs by using Geant 4, and discuss its various implications, taking both of high-luminosity (HL) 
and low-luminosity (LL) GRBs into account. This is because LL GRBs may also make a significant 
contribution to the observed UHECR flux if they form a distinct population. We show that not 
only protons but also heavier nuclei can be accelerated up to ultra-high energies in the internal, 
(external) reverse and forward shock models. We also show that the condition for ultra-high-energy 
heavy nuclei such as iron to survive is almost the same as that for ~ TeV gamma-rays to escape 
from the source and for high-energy neutrinos not to be much produced. The multi-messenger 
astronomy by neutrino and GeV-TeV gamma-ray telescopes such as IceCube and KM3Net, GLAST 
and MAGIC will be important to see whether GRBs can be accelerators of ultra-high-energy heavy 
nuclei. We also demonstrate expected spectra of high-energy neutrinos and gamma rays, and discuss 
their detectabilities. In addition, we discuss implictaions of the GRB-UHECR hypothesis. We point 
out, since the number densities of HL-GRBs and LL-GRBs are quite different, its detemination by 
UHECR observations is also important. 

PACS numbers: 95.85.Ry, 98.70.Rz, 25.20.-x, 14.60.Lm, 96.50.Pw, 98.70.Sa 



I. INTRODUCTION 



Gamma-ray bursts (GRBs) and supernovae (SNe) are 
the most powerful phenomena in the universe. The latter 
is believed to be the origin of high-energy cosmic rays be- 
low the knee w 10 15 ' 5 eV. The former could also acceler- 
ate baryons up to high energies if the dissipation process 
is due to shock dissipation. Waxman, Milgrom & Usov 
and Vietri suggested that ultra-high-energy cosmic rays 
(UHECRs) are produced in GRBs 0, H Eg. Their sug- 
gestions were based on the two arguments. First is the 
requirement that relativistic outflows that make GRBs 
satisfy various conditions for baryons to be accelerated up 
to greater than 10 20 eV. Second is that the energy gener- 
ation rate required to account for the observed UHECR 
flux is comparable to the energy generation rate of ob- 
served gamma rays. The latter argument depends on 
the local GRB rate which is not well known observation- 
ally. If the local cosmological high-luminosity GRB (HL 
GRB) rate is not high enough, which may be suggested 
by recent observations 14, JfJ, [fj , the required baryon load- 
ing becomes larger 0, M,M- ^ n addition, if the UHECR 
spectrum at the source is steeper than that with the spec- 
tral index p ~ 2.0 — 2.3 (as expected in the dip model, 
p ~ 2.4 — 2.7 [{J [HI), the total nonthermal cosmic-ray 
energy of GRBs, necessary for explaining the observed 
UHECR flux, would be larger than the radiation energy. 
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Despite such caveats, the GRB-UHECR hypothesis is 
still one of the most interesting possibilities for expla- 
nation of observed UHECRs, and it should be examined 
by future observations. 

Recent observations have tentatively suggested that 
some GRBs may form a different class (see, e.g., (llTIl^). 
GRBs in this class are called low-luminosity GRBs (LL 
GRBs) or sub-energetic GRBs, which may be more com- 
mon than HL GRBs. If this speculation is true, they can 
provide the significant energy as high-energy cosmic rays, 
neutrinos and gamma rays. Murase et al. 14j suggested 
that cosmic rays can be accelerated up to ultra-high en- 
ergies (> 10 18 - 5 eV) in LL GRBs, and LL GRBs may 
make an important contribution to the observed high- 
energy cosmic-ray flux and the diffuse neutrino back- 
ground [3, EH- Such LL GRBs were associated with 
energetic supernovae called hypernovae. As in cases of 
supernovae, cosmic-ray acceleration is expected at shocks 
formed by the high-velocity ejecta. In fact, if the mag- 
netic field is sufficiently amplified, cosmic-rays can be 
accelerated up to very high energies, even to ultra-high 
energies. Wang et al. [H, [ljj suggested that they could 
explain observed high-energy cosmic rays above the sec- 
ond knee, under the assumption that LL GRBs have the 
enough local rate higher than that of HL GRBs. 

Recently UHECR observations have been greatly ad- 
vanced thanks to (Sourth) Pierre Auger Observatory 
(PAO). The UHECR spectrum deduced by PAO suggests 
the existence of the Greisen-Zatsepin-Kuz'min (GZK) 
cutoff energy which arises from the photomeson produc- 
tion process between UHECRs and cosmic microwave 
background (CMB) photons pjlGjl. In addition, PAO 
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recently revealed the anisotropy and rejected the isotropy 
of UHECRs above 60 EeV at 3 a level [H [S|. It sug- 
gests that the distribution of UHECRs traces the matter 
distribution in the universe, although the origin of UHE- 
CRs is still unknown, that is, AGNs, GRBs and other 
possibilities remain (see the latest review for AGN and 
GRBs, e.g., (22J). Interestingly, the preliminary elonga- 
tion rate data by PAO showed the break around ~ iq 18 ' 35 
eV [23[ . Whether heavier nuclei becomes more important 
above this break as the energy increases (which seems 
conflict with HiRes data) is still under debate. It is be- 
cause there are large uncertainties in hadronic interac- 
tions. On the other hand, recent results on UHECRs 
above 60 EeV by PAO may infer that protons are dom- 
inant sources 120, 

mi, 

although some authors dispute 
[24l [25| . Although the issue of the composition of UHE- 
CRs requires more and more studies, present data seem 
to allow the presence of significant fractions of ultra- high- 
energy (UHE) nuclei around the GZK cutoff energy, even 
if protons are dominant. Therefore, it is one of the im- 
portant issues whether UHE nuclei can survive (i.e., UHE 
nuclei are not photodisintegrated and not depleted due 
to photomeson production during the dynamical time) in 
UHECR production sites. 

In this paper, we study acceleration of high-energy 
cosmic-ray nuclei in very detail by exploiting numerical 
calculations, and discuss various implications for multi- 
messenger astronomy. In Sec. II, we show that cosmic 
rays (protons and heavier nuclei) can be accelerated up 
to ultra-high energies in both of HL GRBs and LL GRBs. 
We consider the internal, (external) reverse and forward 
shock models. 

If cosmic rays can be accelerated up to sufficiently high 
energies, high-energy neutrino production via the pho- 
tomeson production process and/or pp reaction is un- 
avoidable. Now, large neutrino detectors such as IceCube 
KM3Net [23 are being constructed. ANITA [H 
and PAO [Ii| can also detect very high-energy neutrinos. 
We have chances to see such high-energy neutrino signals, 
which arc important as a probe of cosmic-ray accelera- 
tion, in the future. High-energy neutrino signals from 
HL GRBs were predicted in the context of the standard 
scenario of GRBs assuming that observed UHECRs come 
from GRBs [3(1 HH, [H, HU . We also investigated the neu- 
trino emission from GRBs in very detail [3, [3, HH, [H| , 
and various predictions in the Swift era are summarized 
in Ref. [35j . In Sec. Ill, we demonstrate predicted neu- 
trino spectra in some of our models and discuss implica- 
tions for neutrino astronomy. We can see that neutrino 
signals highly suppressed when very heavy nuclei like iron 
can survive. 

In Sec. IV, we also demonstrate gamma-ray spectra 
for some models when heavy nuclei can be accelerated 
up to ultra high energies and survive. In this section, 
we also show that very high-energy gamma rays with 
GeV-TeV energies can escape from the source under con- 
ditions where UHE heavy nuclei such as iron can sur- 
vive. We can expect to detect high-energy gamma-ray 



signals by upcoming Gamma-Ray Large Area Telescope 
(GLAST) and Cherenkov telescopes such as MAGIC as 
long as the source is nearby. Otherwise, high-energy 
gamma rays above ^TcV are attenuated by the cos- 
mic microwave/infrared background (CMB/CIB) pho- 
tons ( 38] and references there in). 

In Appendix F, We also discuss various implications 
of the GRB-UHECR hypothesis (including the HLGRB- 
UHECR, LLGRB-UHECR and Hypernova-UHECR hy- 
potheses) for UHECR astronomy. We see that the deter- 
mination of the source number density of UHECR sources 
and extragalactic magnetic field (EGMF) strength is es- 
sential to test the hypothesis. In this paper, we shall 
hereafter use notations such as Q x = Q/10 x in cgs units. 



II. ACCELERATION OF COSMIC-RAY NUCLEI 
IN GRBS 

In the standard model of GRBs (for reviews, see, e.g., 
[3^ . |40L EH), both of the prompt and afterglow emis- 
sion are attributed to electromagnetic radiation from rel- 
ativistic electrons accelerated in shocks. If relativistic 
outflows that make GRBs contain baryons, not only elec- 
trons but also protons (and heavier nuclei) will be acceler- 
ated. Although the detail of the acceleration mechanism 
is poorly known, we usually assume that the first-order 
Fermi acceleration mechanism works in GRBs, where 
the distribution of nonthermal cosmic rays is given by a 
power-law spectrum under the test-particle approxima- 
tion (see reviews, e.g., HI])- Therefore, we have [H, HU 

p* = e -P, (i) 

where sqr is the energy of cosmic rays in the comoving 
frame, Ucr is the comoving cosmic-ray energy density, 
and p is the spectral index. As a shock becomes highly 
relativistic, the spectral index can deviate from values 
for non-relativistic shocks. In the diffusive small pitch- 
angle scattering regime, p ~ 2.2 is obtained in the ultra- 
relativistic limit [44|. But the lar ge- a ngle scattering can 
lead to harder spectral indices [45l l46l| . For simplicity, we 
shall use p — 2 which is often expected for non-relativistic 
or mildly relativistic shocks with the compression ratio 
r c = 4. The acceleration time scale is written as t acc = 
(v £ N)/(ZeBc), which could be applied to the second- 
order Fermi acceleration mechanism when the Alfven- 
wave speed is close to the light speed [U [H, \^ . Here, 
e jv is the cosmic-ray energy in the comoving frame and rj 
is a pre-factor. For the most efficient acceleration, we can 
expect r\ ~ (1 — 10) [47]. Throughout this paper 77 = I is 
used to obtain the upper limit of the maximum energy. 
Note that 77 ^ 10 would be a more realistic value, and 
the acceleration time scale with r\ = 10 is demonstrated 
in Ref. 0. 

The maximum energy of cosmic-rays is determined by 
several criteria. One is derived from the requirement that 
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the Larmor radius of a particle should be smaller than the 
effective size of the acceleration region (the Hillas con- 
dition). In addition, the acceleration time scale should 
be smaller than the dynamical time scale idyn (which is 
essentially the same as the Hillas condition in our in- 
terested cases) and escape time scale t cac . The maxi- 
mum energy is also limited by the total cooling time scale 
Col = i W7+ i syn+ i ic 1 + t ad 1 - Here, t syn is the synchrotron 
cooling time scale, tic is the inverse-Compton cooling 
time scale and i ac i is the adiabatic cooling time scale. In 
cases of GRBs, we shall use t ad w £ dyn [3(| [H, [5(| [HI- 
t is the photohadronic time scale which includes the 
photodisintegration, photomeson production and pho- 
topair production processes. For protons with sufficiently 
high energies, the photomeson production process is the 
most important, whose energy loss time scale is given by 

c f°° f°° dn 

t vi( e v) = T^I dea pl (e)K p (e)e dee- 2 — , (2) 

where e is the photon energy in the rest frame of pro- 
ton, 7 P is the proton's Lorentz factor, Kip IS the inelas- 
ticity of proton, and £th ~ 145 MeV is the threshold 
photon energy for photomeson production. For heavier 
nuclei than proton, both of the photodisintegration and 
photomeson production processes are important, whose 
time scales are given by the similar expression to Eq. 
© [52j . In order to decide whether a kind of UHE nu- 
clei can survive or not, we calculate the interaction time 
scale of photodisintegration and photomeson pro duction 
by using the numerical simulation kit Geant4 [53[ , which 
includes the cross section data based on experimental 
data [54[. As seen later, it is important to use the accu- 
rate cross section in the high-energy range. Although we 
evaluate t^-y by numerical calculations, the simple ana- 
lytic treatment is often useful. The most frequently used 
approximation for the photomeson production process 
is the A— resonance approximation (see, e.g., [jjQ, HH). 
The corresponding one for the photodisintegration pro- 
cess is the Giant-Dipole-Resonance (GDR) approxima- 
tion. Similarly, we can apply it to the GDR approxima- 
tion for a broken power-law photon spectrum. We have 

where E h N ~ 0.5e r esWArc 2 r 2 /e^ b , a is a photon index 
in lower energies while (3 in higher energies. The pa- 
rameter regions for the upper and lower columns are 
E N < E b N and E N > E b N , respectively. Here e b oh is 
the break energy measured by the observer in the lo- 
cal rest frame, and J7 7 is the total photon energy den- 
sity. For example, a Band function which reproduces 
spectra of the prompt emission, can be approximated 
by a broken power-law spectrum. Spectra of after- 
glows are also expressed by several segments of power- 
law spectra. The cross section at the resonant energy 
e ros is expressed as <7 ros - For GDR of nucleon, we use 



ctgdr ~ 1.45 x 10~ 27 cm 2 A emR ~ 42.65A- 021 MeV 
and Ae ~ 8 MeV [H, HE lHa lH- In the case of 
A— resonance of proton, we use cta ~ 4 x 10 _28 cm 2 , 
e A ~ 0.3 GeV, Ae ~ 0.2 GeV and n p ~ 0.2 [H H||. 
These approximations reproduce numerical results well 
except in high energies. In high energies, effects of non- 
GDR and/or non-A-resonance such as fragmentation and 
multi-pion production can become important, as seen 
later. This tendency can also be seen in t pi . As pointed 
out in Ref. [H, the effect of non-A-resonance such as 
multi-pion production becomes moderately important in 
high energies for spectra of the prompt emission where 
a~ 1. 

For later discussions, let us introduce optical depths 
for photodisintegration of nuclei and photomeson pro- 
duction of protons, and define /at 7 = £dyn/£jV"7 an d 
fpy = idyn/Vy respectively. These quantities express 
whether cosmic rays can survive in the source or not. For 
example, cosmic rays can survive from the photodisinte- 
gration and photomeson production processes if /jv 7 < 1, 
while not if /jv 7 > 1. 

In the following subsections, we shall show that high- 
energy cosmic-ray nuclei can be produced in the internal 
and external shock models for HL and LL GRBs. GRBs 
may be even UHECR sources. In such cases, note that 
cosmic rays have to be accelerated above 10 20 eV from 
observations of the highest energy events and the exis- 
tence of the bump around the GZK cutoff energy. We 
shall allow for possibilities that observed highest UHE- 
CRs arc UHE nuclei. It is because, although the arrival 
distribution and spectrum of UHECRs may be consistent 
with the proton model, some authors claimed that UHE 
nuclei are more important. 



A. Internal Shock Model 

The internal shock model is one of the most frequently 
discussed models in order to explain the prompt emission. 
Relativistic outflows make internal collisions, which lead 
to internal dissipation via shocks. The formed shocks will 
be mildly relativistic shocks, and charged particles will 
be accelerated at those collisionless shocks by some mech- 
anism such as the Fermi acceleration mechanism. HL 
GRBs have luminosity L 1 ~ 10 51-52 ergs/s and the ob- 
served peak energy ~ 10 2,5 keV. The isotropic equiv- 
alent energy is E™° ~ 10 53 ergsZ 7)52 [£t-i/(l + z)]N 2 . 
Here 5t is the variability time which can vary in the 
very wide range and N is the number of collisions. A 
typical collision radius will be expressed by commonly 
used relation, r m 1.2 X 10 14 r| 5 [<5t/0.02(l + z) s] cm. 
Of course, this radius has to be smaller than the de- 
celeration radius where the afterglow begins. As dis- 
cussed in Refs. [3(| [58|, an internal collision radius 
is one of the most important quantities for the pho- 
tomeson production. A collision radius r ~ 10 13 ~ 15 
cm is frequently assumed. The magnetic field is given 

by b = 7.3 x 10 4 G4 /2 (r sh (r sh - 1) / 2 ) 1/2 x^r^-/ . 
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FIG. 1: The acceleration time scale and various cooling time 
scales of proton and oxygen in the internal shock model for HL 
GRBs. Energy and time scales are measured in the comoving 
frame of the outflow. Used parameters are Lt = 10 5 
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10 2 ' 5 keV, T = 10 25 , r = 10 14 cm and £ s (f« e s /e e 



Thick lines show numerical results on the photomeson and/or 
photodisintegration time scales. Thin lines show analytic re- 
sults obtained by the resonance approximation. In the high 
energies, the effect of the non-resonant cross section becomes 
important. Note that this parameter set implies that a sig- 
nificant fraction of the energy carried by protons goes into 
neutrinos. 



Here, Lm is the kinetic luminosity of outflows and r sn is 
the relative Lorentz factor between two subshells. The 
typical width of subshells in the comoving frame I is typ- 
ically given by I = r/T. 

Next, let us evaluate maximum energies of cosmic rays 
by using several criteria. First, t acc — tdyn ~ tad leads to 
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Note that the Hillas condition = I = r/T is already 
satisfied. Second, t acc 



t syn leads to 
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Therefore, we can expect that cosmic rays can be acceler- 
ated up to ultra-high energies in the internal shock model 
of HL GRBs unless other cooling time scales such as t P7 
are important. 

The inverse Compton cooling time scale can be also 
calculated. For evaluation of tic, we need to give a 
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FIG. 2: The same as Fig. I. But used parameters are Lt = 
10 51 ergs-\ e b oh = 10 2 ' 5 keV, F = I0 3 , r = 10 15 cm and 
Cb(~ es/e e ) = 1. Thick lines show numerical results with 
the synchrotron self-absorption cutoff 1 eV (in the comoving 
frame). Thin lines show analytic results without the low- 
energy cutoff for comparison. We can see that the effect of 
self-absorption is not so important for the maximum energy 
due to the effect of the non-resonant cross section in the high 
energies. Note that this parameter set allows for survival of 
UHE nuclei. 



photon spectrum. The photon spectrum for the prompt 
emission is often approximated by the broken power-law 
as 
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where is the luminosity at the break energy measured 
by the observer in the local rest frame, e min is the mini- 
mum cutoff due to synchrotron self-absorption and e max 
is the maximum cutoff due to pair-creation in the co- 
moving frame. In this paper, we set e mva = 1 eV and 
£ max _ 20 MeV. In this paper, we shall use a = 1 and 
[3 = 2.2 as photon indices. Although tic is calculated, 
we can usually ignore this cooling time scale due to the 
Klcin-Nishina suppression. Hence, E™* c is usually larger 
than E^ n . 

The effect of the photomeson production process of 
protons is investigated in detail in Ref. [if]. For details, 
see Refs. [1, [35j] and references there in. At smaller colli- 
sion radii, t^> 7 can be more important than other cooling 
time scales such as t syn . UHE protons are not depleted 
only at sufficiently large radii. The cross section of pho- 
todisintegration is larger than that of photomeson pro- 
duction, so that survival of UHE nuclei is more difficult 
than survival of UHE protons. We evaluate the maxi- 
mum energy due to photomeson production and/or pho- 
todisintegration, E™* and/or -E™ a ^ 7 , from t acc = t p7 
and/or t acc = t Nj . 

Our numerical results on various time scales in the in- 
ternal shock model for HL GRBs are shown in Figs. 1 
and 2. For calculations, we give the total photon energy 
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/ dee(^). The magnetic field is 
{e B /e e )U T 

.AT - 10 2..' 



density U~ ( by f/ 7 

given by Ug = S,bU^ » (es/e e )LL. Fig. 1 shows the 
result for r = 10 14 cm and T = 10 . Thick lines show 
numerical results while thin lines show curves when we 
use resonance approximations. At sufficiently high ener- 
gies, we can see that effects of the cross section in the 
non-resonance region become important. For this typical 
parameter set, we have the effective optical depth for the 
photomeson process f pl ^0.3. Hence, the efficient neu- 
trino production occurs in this parameter set, which can 
lead to detectable neutrino signals d,|3(j. However, we 
cannot expect survival of UHE nuclei in such cases. Al- 
though they can be accelerated up to very high energies, 
UHE nuclei cannot survive in the sense that /n^ > 1. 

Fig. 2 shows the result for r = 10 15 cm and T = 10 3 . 
In this parameter set, we expect that UHECR nuclei can 
survive. Protons and oxygens can be accelerated up to 
~ 10 20 eV and ~ 10 21 eV, respectively. This is just be- 
cause the photon density becomes small enough at large 
radii. Note that we have the effective optical depth for 
the photomeson process / P7 ~ 10~ 3 , which is a very small 
value. In this parameter set, although we can expect sur- 
vival of UHE nuclei, the magnetic field is also expected to 
be rather weak as long as we use the conventional value 
£b = 1- Here, we have B ~ 5.8 x 10 2 G, which seems to 
be somewhat insufficient to explain the observed break 
energy by the optically thin synchrotron model. (Then, 
we expect e^ b ~ 12keVr 2 . 5 (r sh - lfe 2 e g 2 f~ 2 B 3 , where 
g = g(p) — (p — 2)/(p— 1) and p is the electron spec- 
tral index). But the detailed mechanism of the prompt 
emission has not been revealed yet 59], so that we do 
not consider this problem in this paper. In Fig. 2, we 
also show curves for spectra with the minimum cutoff en- 
ergy (thick lines) and that without the cutoff (thin lines). 
Wang et al. [T3| suggested that the effect of minimum 
cutoff energy (due to self-absorption) can increase the 
maximum energy of heavy nuclei. However, we see that 
the effect of self-absorption is not so important even at 
the highest energies due to the effect of the cross section 
in the non-resonance region, especially where the frag- 
mentation process is important. 

Next, let us consider LL GRBs such as GRB 060218 
and GRB 980425. For example, GRB 060218 has low 
luminosity L 1 ~ iq46-47 er g S /s, which is much smaller 
than that of usual HL-GRBs [60] ■ The duration time 
is T ~ 3000 s, and the isotropic equivalent energy is 
E™° ~ 10 50 ergs. The observed peak energy is e b oh ~ 5 
keV, hence this event is classified as a x-ray flash. Various 
interpretations for GRB 060218 exist. Several authors 
suggested the shock break-out model [fJl], |6^|. Other 
authors argued that this event can be explained by the 
internal-external shock model [H, HU, [65| . In this paper, 
we adopt the conventional internal shock model follow- 
ing Toma et al. and we shall adopt the Lorentz factor 
r ~ (5 — 10), which is suggested in Refs. [H, [65[. 

Fig. 3 shows one of our results for r = 6x 10 15 cm and 
r = 10, which correspond to 8t ~ 1000 s when we use 
r w 2T 2 cSt. Results of other parameter sets are presented 
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FIG. 3: The acceleration time scale and various cooling time 
scales of proton, helium, oxygen and iron in the internal shock 
model for LL GRBs. Energy and time scales are measured 
in the comoving frame of the outflow. Used parameters are 



L b = 1.5 x 10 4b ergs- 



pi 

e ob 



= 5 keV, r = 10, r = 6 x 10 1 



cm (corresponding to T — 10 s) and £s(~ es/e e ) = 1. Note 
that this parameter set allows for survival of UHE nuclei. 
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FIG. 4: The maximum energy of cosmic rays (thick lines) and 
the optical depth for photomeson production or photodisinte- 
gration (thin lines) as a function of the collision radius. Used 
parameters are Lb = 1.5 x 10 46 ergs -1 , e b oh — 5 keV, T = 10 
and £b(~ es/e e ) = 1 for LL GRBs. The comoving shell width 
is set to I = r/r. Note that UHE nuclei can be produced in 
LL GRBs for r > 10 15 cm. There, UHE nuclei can survive. 
At sufficiently large radii, _E^ ax becomes -B^ ax d . 



in Appendix A. In this parameter set demonstrated in 
Fig. 3, we can expect that not only UHECRs can be 
produced but also accelerated UHE nuclei can survive. 
This is just because the photon density is small enough 
at large radii. Note that we have the effective optical 
depth for the photomeson process f pl ~2x 10~ 3 , which 
is a small value. In this parameter set, although we can 
expect survival of UHE nuclei, the magnetic field seems 
to be somewhat weak when we use the conventional value 
£b = 1. In Fig. 3, Fe nuclei can be accelerated up to 
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~ 10 21 eV, while protons only up to ~ 10 19 5 eV. 

As noted before, one of the important parameters is 
a collision radius. Hence, we show maximum energies of 
proton and heavy nuclei £/^ ax as a function of a collision 
radius r in Fig. 4. In the same figure, f pl and /at 7 
are shown. We can see that the maximum energy of 
proton is determined by the photomeson production at 
sufficiently small radii, but usually determined by the 
synchrotron cooling and adiabatic cooling processes at 
large radii (r > 10 12 ' 5 cm). The maximum energy of 
nucleon is usually determined by the photodisintegration 
and adiabatic cooling. Survival of UHE nuclei is possible 
only at r > 10 15 cm in this parameter set. 



B. Reverse Shock Model 

The expanding fireball strikes the surrounding medium 
and will form a reverse shock and forward shock. The 
shocked ambient and ejecta materials are in pressure bal- 
ance and are separated by a contact discontinuity. In the 
original standard model, the reverse shock is thought to 
be short-lived, which exists during the initial decelera- 
tion of the fireball. During this phase, optical/infrared 
flashes were expected. Indeed, some optical flashes can 
be interpreted as the reverse shock emission. However, 
recent observations have reported the tentative lack of 
bright optical/infrared flashes [66| . Although it is one 
of the open problems in the Swift era, we do not man- 
age this problem in this paper. Recently, the modern 
version of the reverse shock model was developed to ex- 
plain the shallow decay emission in the early afterglow 
phase [13, l68| . In such models, the plateau shape can 
be achieved by requiring the appropriate distribution of 
Lorentz factors of the ejecta and the suppression of the 
forward shock emission. For the reverse shock emission 
to emerge in the x-ray band, we can consider a number 
fraction of the shocked electrons that are injected into 
the acceleration process /J to be smaller than the unity. 
In this paper, we shall consider cosmic-ray acceleration in 
the early afterglow phase under the reverse shock model 
developed for the shallow decay emission, and adopt the 
small f T e following Refs. |E|6|. 

The cosmic-ray production in the reverse shock model 
was discussed in Refs. [!, |31 L 
study can be found in Ref. 



32] 



35] 



The more detailed 
As demonstrated in 
these references, cosmic-ray acceleration up to ultra-high 
energies at a reverse shock of HL GRBs is possible. In 
this paper, let us demonstrate that UHE cosmic-ray pro- 
duction at a reverse shock of LL GRBs is also possible. 
First, t acc = tdyn ~ tad leads to 



(1 + z)ETZ d 



ZeB r A r x 



4.7 x 10 20 eVZr)' 1 



x 4 2 El!lA'*Tl'\ (7) 



where T is the duration time measured by the observer 
in the local rest frame. Here, we have assumed the thick 



ejecta case, where the total ejecta thickness ~ cT is larger 
than the shocked region at the crossing time ~ r x /2TQ. 
Here Tq is the initial Lorentz factor. In the thin ejecta 
case, the corresponding expression can be derived easily. 



Second, i acc = t syn leads to 



(1 + z)E_ 



max 
7V,syn 



QitZe 



N 



V Z i axB , x ri m e 
~ 1.3 x 10 21 eVA 2 Z- 3 /V 1/2 

X e B Ay, 50 n 2 1 3 ■ ( 8 ) 

Therefore, we can expect that UHECR production at a 
reverse shock is possible not only for HL GRBs but also 
LL GRBs. Note that even protons could be accelerated 
above 10 20 eV unless es is too small. 

Photon spectra in the reverse shock model can be cal- 
culated by exploiting the reverse-forward shock model. 
For example, in the slow cooling regime (e m < e sa < e c ), 
we obtain 



dn 
de 



where [351 Kf 



-(e/e 



(e sa /e m )~— {£ m /£ sa )^{e/e m y 

(e sa / £ m )~ £ * i ( £ /e sa )5 
p+i 

(e/£ m y— 

p+1 p+2 

{e c /e m ) ~{e/e c )~~ 

(9) 



(10) 



Anri ce r ' 



Here e n = mm[£ m ,e c ], and e m , e c and e sa are the injec- 
tion, cooling and self-absorption energy in the comoving 
frame, respectively. L E ^ raax is the comoving specific lumi- 
nosity per unit energy at the injection or cooling energy. 
From Eq. ©, we can calculate t\c and tjv 7 - 

We show the numerical result of the reverse shock 
model for LL GRBs in Fig. 5. The used parameters 
are E c] = 2 x 10 51 ergs, n = 10 2 c m - 3 , e r e = e r B = 0.1, 
fl = 0.025, T = 3000 s and the initial Lorentz factor 
r = 5. As seen in Fig. 5, the UHECR production 
is possible, where protons and Fe nuclei can be acceler- 
ated up to ~ 10 20 ' 2 eV. But we cannot expect survival 
of UHE nuclei in this parameter set due to the copious 
photon field. For UHE nuclei accelerated above 10 20 eV 
to survive, for example, much smaller e r e is needed. If we 
adopt such parameters, we can expect both UHE protons 
and nuclei in this model. 



C. Forward Shock Model 

The expanding fireball strikes the surrounding medium 
and will form a forward shock. The self-similar behav- 
ior after the deceleration time t& ec is described by the 
Blandford-Mckee solution. The standard external shock 
model based on this solution can reproduce observations 
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4 6 

log(s N [GeV]) 

measured in the comoving frame 



FIG. 5: The acceleration time scale and various cooling time 
scales of proton, helium, oxygen and iron in the reverse shock 
model for LL GRBs. Energy and time scales are measured 
in the comoving frame of the outflow. Used parameters are 
E ci = 2 x 10 51 ergs, n = 10 2 cm -3 (ISM), r = 5, T = 3 x 10 3 
s, e r B = 0.1, e r e = 0.1, fl = 0.025 and p = 2.4. Note that this 
parameter set does not allow for survival of UHE nuclei. 




log(t [s]) 



at the late time. Although the original forward shock 
model actually fails to explain the shallow decay emis- 
sion^ 9 we consider the original forward shock model for 
simplicity. By exploiting the theory, we immediately have 
r oc t~«, r oc t~£ and B oc £~s. By using them, we can 
derive the time-dependence of various quantities. Typ- 
ical parameters for HL GRBs are ~ i0 52 " 53 ergs, 

2. In 



lcm- 3 (ISM), e B ~ 0.01, e e 



0.1 and p 




5 

log(t [s]) 



FIG. 7: The same as Fig. 6. But the wind-like circumburst 
medium is assumed (A* = 0.77). UHE nuclei can be pro- 
duced in a forward shock of LL GRBs. Survival of UHE nu- 
clei becomes easier with time. At sufficiently late time, iJ™ ax 



obtaining above parameters, all the electrons are often 
assumed to be accelerated. However, this may not be 
true. As Eichler & Waxman [69[ discussed, only a frac- 
tion f e of electrons may be accelerated, that is, f e can 
be smaller than the unity. When we consider the small 
/ e , the kinetic energy of ejecta can be larger by l// e . 
Such small values of f e will be tested by observations of 
polarization in the future [toI ] . 

Cosmic-ray acceleration and neutrino production in 
the forward shock model for HL GRBs were discussed 
in detail by Dermer [H, |50T| . Hereafter, let us demon- 
strate UHECR production at a forward shock for pa- 
rameters of LL GRBs. In this paper, we shall adopt 
2 x 10 51 ergs, n = 10 2 cm- 3 (ISM), e B = 0.01, 
= 0.01, f e = 0.1 and p — 2. This parameter set is 



FIG. 6: The maximum energy of cosmic rays (thick lines) and 
the optical depth for photomeson production or photodisin- 
tegration at 10 20 eV (thin lines) as a function of the time tion t a 
after the burst. Used parameters are E e j — 2 x 10 51 ergs, 
n = 10 2 cm -3 (ISM), e B = 0.01, e e = 0.01, f e = 0.1 and 
p — 2. The jet-break time is set to tj = 2 x 10 4 s. Note that 
the deceleration time is t& cc ~ 8000 s. UHE nuclei can be 
produced in a forward shock of LL GRBs. Survival of UHE 
nuclei becomes the most difficult at the jet-break time. At 
sufficiently late time, i5^ ax becomes E™* d . 



consistent with parameters obtained by Toma et al. [65( 
except that we use f e = 0.1. 

Next, let us estimate maximum energies of cosmic rays. 
First, the maximum energ y d etermined from the condi- 



tdy 



t. 



ad 



48. 4£ 



(1 



\ pmax 



ZeTB 



ISM. 



7.1 x lO^eVZ??- 1 



B 



ISM 



2 3/ V 



(11) 



Here (1 + z)t is the time after the prompt emission 
measured by the observer. When we use the conven- 
tional value of the ISM magnetic field, we cannot ex- 
pect UHECR production at a forward shock of GRBs. 
However, the upstream magnetic field might be ampli- 
fied significantly via some process, e.g .. the non-resonant 
cosmic-ray streaming instability [33|, [L!j (but see also 
Ref. |z3). Instead, if the second-order Fermi accelera- 
tion works well [7|[ , we could use the magnetic field in the 
downstream, which is likely to be amplified, as indicated 



8 



by observations. If the upstream and/or downstream 
magnetic fields are amplified up to near the equiparti- 
tion value, i acc = id yn ~ iad leads to 

(I + z)E^ d = ZeBfrr,- 1 ~ 2.7 x 10 20 eV Z^ 1 
The maximum energy limited by synchrotron cooling is 



(I + z)E$%L 



6irZe Tm 2 N c 2 



~ 2.5 x 10 20 eVA 2 Z- 3 /V 1/2 

x ^Elif^h^. (13) 



From Eqs. flT2J) and fig]), we can expect that UHECR 
production is also possible at a forward shock of LL 
GRBs. 

A photon spectrum in the forward shock model can be 
also calculated from the theory. For example, in the slow 
cooling case (e sa < e m < e c ), we have 

(e sa /e m Y^{e/e sa f 
dn _ J L e ,maxe- (e/emax ) J (e/e" 1 )" 1 
de 



47rr 2 c£ I 



p+i 

2 

P+l 



(e/e T 

(e7£ m )"^(e/e c ) _iT " 



(14) 



Here, = Te n 

nl7 „ V „2 f -2 1/2 9 pl/2 ,-3/2 c 
U - 17eV 3-1 Je,-l e S,-2 e e -1^,51*4 i £ ob 
1 n i7 -3/2 t-,-1/2 -1,-1/2 , T 

l.OeV e B L 2 E ci ^ n 2 i 4 and L eob , max = 

7.8 x lO 56 s- 1 /e,-i4 /2 -2^cj,5in2 /2 (0 P /O.6), where </> p is 
an order-of-unity factor calculated by Wijers & Galama 

The photodisintegration and/or photomeson time 
scales are evaluated from Eq. @. Here, let us estimate 
these time scales analytically in the slow cooling case by 
applying Eq. ([3]). After replacing e b with s c , resonance 
approximations give 



t N \ ~ g^ 1 { ( t >p /Q.&)e e ^lXlli^ /S h b/ \E N /E b N ) 



7b \V2 



1.1 x 10 8 s 1 (for A-resonance) 

5.7 x lCT 8 /! 121 s" 1 (for GDR) 



(15) 



The above expressions are obtained for En > E h N = 
0.5£ r esWi./vc 2 r 2 /££ b . Hence, we have /jv 7 oc t» and 
t~i for ^jv > E b N . Of course, the time- 
dependence differs for En < E h N . In the ISM case, we 
can see that //v 7 increases with time. Note that the tem- 
poral index of Jn-^ is different from that used in Wang 
et al. [13], because we consider En ~ 10 20 eV, where 
En > E b N is typically expected rather than En < E b N 
considered by them. 

In the wind- medium case, we have E™ 3 ^ x t~ J an< ^ 



E 



N ,syn 



is. Furthermore, we can derive E 



N,N-y 



OC 



and /at 7 oc i -1 for En > E b N . In the wind- medium case, 
we can see that /jv 7 decreases with time. It is because 
the ambient density drops as r~ 2 , which leads to rapid 
drop of the photon density with time. 

In the light curve of afterglows of HL GRBs, we can 
often find the jet break around tj ~ 10 5 s. After the jet 
break, we have r oc i°, T oc i -1 / 2 and B oc i -1 / 2 . There- 
fore, for example, we obtain E^ a ^ A oc £~3 ( Effi^ oc £~3, 

e n*n~i « i° and /at 7 oc i~5 for > £^ in the ISM 
case. After the jet break, the photon density will become 
smaller and smaller. Hence, we expect photodisintegra- 
tion and photomeson production become less important. 
Similarly, we can obtain the time-dependence of various 
quantities for the wind-medium case. In both cases, the 
effect of the jet break makes survival of UHE nuclei eas- 
ier. It has not been pointed out in previous studies. 

Results in the forward shock model for LL GRBs are 
shown in Figs. 6 and 7. In Fig. 6, we show maximum 
energies of cosmic rays E^ ax and the optical depth for 
photodisintegration /jv 7 at En = 10 20 eV as a function 
of time i. From this figure, we can expect UHE nuclei can 
be produced at a forward shock of LL GRBs. For heavy 
nuclei (O and Fe), the photodisintegration determines the 
maximum energy at the earlier time, but the adiabatic 
cooling becomes more important at the later time. For 
light nuclei, the adiabatic cooling is almost always the 
most important. Note that the jet-break effect can be 
important for survival of UHE nuclei, as expected. In 
fact, after the jet break, /p e7 can decrease with time and 
become smaller than the unity. 

In Fig. 7 we show the results for the wind-medium 
case, In this case, we also expect acceleration and sur- 
vival of UHE nuclei. Since the photon density decreases 
with time, the maximum energy of heavy nuclei, which 
is determined by E^ 3 ^ , can increase until the adiabatic 
cooling becomes more important. As a result, UHE nu- 
clei could be accelerated up to > 10 20 eV in this param- 
eter set. 



III. IMPLICATIONS FOR NEUTRINO 
ASTRONOMY 

Sufficiently high-energy cosmic rays cannot avoid the 
photomeson production process, where the pion produc- 
tion threshold energy is 145 MeV. Generated pions, kaons 
and other mesons can produce high-energy neutrinos. 
Such high-energy neutrino signals, if detected, are very 
important as a direct probe of cosmic-ray acceleration. 
High-energy neutrino emission from accelerated protons 
at shocks was predicted by Waxman & Bahcall [30( in the 
internal shock model. Predictions were also done in the 
reverse- forward shock model [3ll [32|, [33j . Further predic- 
tions in the Swift-era are found in Refs. [lH l34l |35| . 

The important quantity for the neutrino flux originat- 
ing from protons is the effective optical depth for the pho- 
tomeson production process / P7 , which also represents an 
energy fraction of protons carried by mesons, as long as 
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f P1 < 1. When fjyy > 1, accelerated protons are depleted 
by photomeson production, which also means the effi- 
cient photomeson production. Such conditions are sat- 
isfied when internal shocks making the prompt emission 
occur at sufficiently small radii [8]. In addition, / P7 > 1 
is also possible in the case of flares and late prompt emis- 
sions [H, HH. 

On the other hand, when / P7 < 1, high-energy pro- 
tons can survive without complete depletion. Such con- 
ditions are satisfied when internal shocks making prompt 
emission occur at sufficiently large radii. In the reverse- 
forward shock model (with the ISM environment), / P7 < 
1 is usually expected except in the highest-energies. How- 
ever, we still expect that a significant fraction of the non- 
thermal proton energy is released as neutrinos unless / P7 
is too small [U HH, Ha] • 

In most previous works, all the cosmic rays are as- 
sumed to be protons. However, heavier nuclei may be 
entrained, which can also be accelerated up to ultra-high 
energies as shown in Sec. II. We can relate f pl with /jy 7 
as 



f P -f(E N ) 
{f Nl {E N )/A) 



Q.2A- 




NJ 



(16) 



Here, the parameter regions for the upper and lower 
columns are {E h N <) Em < E h v and En > E b , respec- 
tively. From Eq. (fT6|) . we can see that / P7 is small at 
radii where UHE heavy nuclei can survive, /at 7 < 1. 

We can calculate neutrino spectra by using the same 
method as in Ref. Q. When UHE heavy nuclei such 
as iron can survive, we easily expect that neutrino emis- 
sion is not efficient from Eq. (fTS)) . For calculations, we 
assume that cosmic rays are all protons in cases where 
UHE irons cannot survive, while proton 75 % and iron 
25 % in cases where UHE irons can survive. Although 
this cosmic-ray composition is just ad hoc, our results on 
neutrino spectra in GRBs do not depend on the composi- 
tion so much, especially when UHE heavy nuclei survive, 
fNf < I- The reason is as follows. The cross section of 
photomeson production for heavy nuclei is roughly ex- 
pressed as cr m eson — Acta (rather than A 2 / 3 cta which is 
given in Ref. lTf), while the inelasticity around the A- 
resonance is around 0.2/ A. Hence, the energy loss time 
scales due to the photomeson production of heavy nu- 
clei is written as t mcson i P7 . Hence, protons and irons 
generate similar neutrino fluence levels in cases of GRBs, 
as long as both are accelerated up to sufficiently high 
energies [zS[76j|. 

As an example, let us show neutrino spectra expected 
in the prompt emission phase, which is the most fre- 
quently discussed since the prediction by Waxman & 
Bahcall [3(J. We can evaluate f pl and / mcson analyti- 
cally by using the A-resonance approximation (30l . l34j . 
For HL GRBs, we obtain 



/p7 ~ fn 



0.3 



J b,51.5 



a — 1 



ri 4 r 2 . 5 2 4 b! 3 16kcV [ (E p /E b p ) 



Here, the parameter regions for the upper and lower 

„, respectively, and 
2.5 due to the effect of 



columns are E p < E h and E p > E b 



we have multiplied a factor of ' 
the multi-pion production which is important for a ~ 1 
spectra 1341. For LL GRBs, we can apply the above ex- 
pression [14J and have 



fp~( ~ fn 



1.4x10" 



Jfc,46.2 



^i5. 8 ri 2 ^ b 



5kcV 



(E P /E b p ) 
(E P /E b p ) 



8-1 



(17) 



Here, the parameter regions for the upper and lower 
columns are E p < E b and E p > E b , respectively. Our 
results are shown in Fig. 8. In fact, the above analytic es- 
timations agree with numerical results. For example, let 
us consider parameter sets demonstrated in Fig. 1 for HL 
GRBs and Fig. 3 for LL GRBs. For the former set with 
the source redshift z = 0.1 (E™° — 10 53 ergs and £ acc = 
20), we have E 2 v (j) v ~ {l/A)f P7 E 2 {dN p so / dE p ) / {AitD 2 ) ~ 
3 x 10~ 4 ergcm -2 , which agrees with the thick solid line 
shown in Fig. 8. For the latter set with the source 
redshift z = 0.005 (E™° = 10 50 ergs and £ acc = 10), 
we have E%<f> v ~ (l/4)f P7 E 2 (dN p BO /dE p )/(4nD 2 ) ~7x 
10~ 7 ergcm -2 , which also agrees with the thin dashed 
line shown in Fig. 8. Note that such low redshift bursts 
(at ~ 20 Mpc) have not been observed yet (e.g., ~ 40 
Mpc for GRB 980425). But we may see such bursts if LL 
GRBs occur in e.g., Virgo cluster. The expected muon 
event rates by IceCube are also shown in the figure cap- 
tion of Fig. 8. As stressed in the previous paragraph, 
survival of UHE heavy nuclei means that neutrino emis- 
sion is inefficient, so that it would be difficult to expect 
detection of neutrino signals by near-future neutrino tele- 
scopes such as IceCube. 

Since it is difficult to see neutrino signals from one 
GRB event, we may need to see many neutrino events 
as the cumulative neutrino background. As we can see 
from Eqs. (|C1[) and (|C2[) . the cumulative neutrino flux 
can be estimated from min[l,/ P7 ] and a given cosmol- 
ogy (see Appendix C). We typically expect min[l, / P7 ] ~ 
(0.01 — 1), for example, in the internal shock model 
for HL GRBs with T < 10 2 5 and r < 10 15 ' 5 cm. 
Smaller values are possible only at larger radii and/or 
for larger Lorentz factors. Survival of UHE heavy nu- 
clei such as iron requires such relatively extreme pa- 
rameter sets, which leads to f pl ~ 10~ 3 . As a result, 
the expected cumulative neutrino flux under the GRB- 
UHECR hypothesis is E 2 <§ v - 10" 8 GeVcm^s^sr- 1 
for the parameter set demonstrated in Fig. 1, while 
E 2 ^^ ~ 3xl0~ n GeVcm _2 s _1 sr _1 for the parameter set 
demonstrated in Fig. 2. The corresponding muon event 
rates by IceCube are N p ~ 50 events/yr and ~ 0.05 
events/yr, respectively. Since the neutrino flux from nu- 
clei is very similar to that from protons when acceler- 
ated heavy nuclei survive, we can use results obtained in 
Murase & Nagataki for mixed composition cases where 
UHE nuclei can survive. The detailed numerical calcula- 
tions on the cumulative neutrino background are found 
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FIG. 8: Energy fluences of neutrinos from one nearby GRB 
event. Solid lines and dashed lines show HL GRB with _B' so = 
10 53 ergs at z = 0.1 and LL GRB with E™ = 10 50 ergs at 
z = 0.005, respectively. A thick solid line shows the HL- 
GRB neutrino spectrum for r = 10 14 cm and F = 10 2 ' 5 where 
heavy nuclei cannot survive, while a thin solid line shows the 
HL-GRB neutrino spectrum for r = 10 15 cm and F — 10 s 
where heavy nuclei can survive (see Figs. 1 and 2). A thick 
dashed line shows the LL-GRB neutrino spectrum for r = 

9 x 10 14 cm and T = 10 where heavy nuclei cannot survive, 
while a thin dashed line shows the LL-GRB neutrino spectrum 
for r — 6 x 10 15 cm and F = 10 where heavy nuclei can 
survive (see Figs. 3 and 10). The cosmic-ray composition 
with proton 100 % is assumed for thick lines, while proton 
75 % and iron 25 % for thin lines. The nonthermal baryon 
loading factors £ acc = Ucn/U-, are set to 20 for HL GRBs and 

10 for LL GRBs, respectively (see Appendix B). We also use 
£b = Ub /Ury = 1. Expected muon event rates by IceCube are 

~ 1 event for the thick solid line, 7V M ~ 0.001 events for 
the thin solid line, iV M ~ 0.2 events for the thick dashed line 
and ~ 0.002 events for the thin dashed line. 

in Refs. [1, [13, [H], HU . In Ref. [1] , neutrino spectra are 
shown for various collision radii and it is useful to com- 
pare set A and set B in Figs. 15-17, for example. So far 
we have considered the internal shock model. For other 
models, see Appendixes D and E. 



IV. IMPLICATIONS FOR GAMMA-RAY 
ASTRONOMY 

Not only neutrinos but also high-energy gamma rays 
originating from cosmic rays (cosmic-ray synchrotron 
radiation), neutral pions, and muons, electrons and 
positrons from charged pions will be produced. However, 
such high-energy gamma rays generally suffer from the 
internal attenuation processes, especially in the internal 
shock model, as discussed in many papers (see, e.g., [36[ 
and references there in). The copious photon field also 
plays an important role on the efficient photomeson pro- 
duction, so that we cannot expect that GRBs are bright 
in ^TeV gamma rays when bright in neutrinos (see Refs. 



[371 l77j and references there in). In other words, when 
fpy becomes small enough, we can expect that the opti- 
cal depth for pair-creation / 77 becomes smaller than the 
unity (hence high-energy gamma rays from far sources 
will be attenuated by CMB/CIB photons rather than 
seed photons in the source) . As seen in the previous sec- 
tion, when UHE heavy nuclei can survive, / P7 is much 
smaller than the unity. Therefore, we may expect escape 
of very high-energy gamma rays from the source when 
UHE heavy nuclei can survive. Such high-energy gamma- 
ray signals are useful as signatures of UHECR accelera- 
tion, so that they are important although the distinction 
between hadronic and leptonic components is not so easy 
in general. Note that, when / P7 is large enough, one may 
not expect high-energy gamma rays due to internal atten- 
uation. But proton-induced cascaded gamma-ray signals 
could be seen, e.g., when the baryon loading factor is 
large enough £ acc ^> 10 [z3, [78[. Hence, observations of 
high-energy gamma rays arc still important as a probe of 
the UHECR acceleration. See also |79j |. 

As an example, let us consider the internal shock model 
with a = 1 and (3 — 2.2. The optical depth for pair- 
creation process (which is usually relevant for escape of 
high-energy gamma rays) can be app roximately written 
as / 77 (£ b) — 0-12<7T£(dn/de)~l [UHO]. Here e is de- 
fined as (m e c 2 ) je. When a ~ 1, the optical depth for 
pair-creation / 77 roughly reaches the maximum value at 

~ e b ob = T 2 (m e c 2 f /s b ob . (Note that, in fact, / 77 (e Q b) in- 
creases logarithmically above ~ e b oh while decreases above 
e„b-) O n t ne other hand, as seen in the previous section, 
/jV7 is also determined by the photon field. Therefore, 
we can relate the two quantities. Note that, in cases of 
a ~ 1, the numerically calculated //v 7 will be effectively 
larger than that evaluated analytically by using Eq. Q 
by a factor of ~ 3. Taking into account this factor, / 77 
at e Q b can be written as 

Here, the parameter regions for the upper and lower 
columns are e b < s b oh and e b > £ob' respectively. From 
the above expressions, we can expect that the cosmic- 
ray source where UHE irons can survive (/Fe-y ^ 1 ) is 
almost completely thin for pair-creation (/ 77 < 1). This 
implies that very high-energy gamma rays are expected 
in the prompt phase, if UHE irons can be accelerated 
at internal shocks and survive. In Fig. 9, we demon- 
strate gamma-ray spectra from one HL GRB at z = 0.1 
and LL GRB at z = 0.005. Gamma-ray spectra given 
by Eq. (6) and originating from cosmic-ray synchrotron 
radiation are shown. The former usually attributes to 
synchrotron (or jitter) radiation from accelerated elec- 
trons. The latter component should exist if cosmic rays 
can be accelerated to ultra high energies. Very high- 
energy gamma rays can be also produced via synchrotron 
self-inverse Compton scattering, which should contami- 
nate the cosmic-ray synchrotron component. Whether 
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FIG. 9: Energy fluences of gamma rays from one nearby 
GRB event. Dashed and dotted lines show HL GRB with 
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ergs at z = 0.1 and LL GRB with E 1 ™ = 10 JU 
ergs at z — 0.005, respectively. A thick dashed line shows the 
HL-GRB gamma-ray spectrum for r = 10 15 cm and V — 10 s 
where heavy nuclei can survive (see Fig. 2). A thin dashed 
line shows the only cosmic-ray synchrotron component from 
one HL GRB without attenuation by the CMB and CIB. A 
thick dotted line shows the LL-GRB gamma-ray spectrum for 
r — 6 x 10 15 cm and T = 10 where heavy nuclei can survive 
(see Fig. 3). A thin dotted line shows the only cosmic-ray 
synchrotron component from one LL GRB without attenua- 
tion by the CMB and CIB. Note that synchrotron self-inverse 
Compton components by accelerated electrons are not shown. 
We assume the comsic-ray composition with proton 75 % and 
iron 25 %. £ aC c = Ucr/Uj and £b = Ub/U~, are set to 20 
and 10 for the HL GRB case, while 10 and 1 for the LL GRB 
case (see Appendix B). We use the low-IR model for the CIB, 
which is presented by Kneiske et al. [8lj . Fluence sensitivity 
curves of GLAST and MAGIC are also shown SI III. 



such a synchrotron self-inverse Compton component ex- 
ists or not, we can expect to detect high-energy gamma- 
rays by GLAST and/or MAGIC in the future, as long as 
accelerated UHE heavy nuclei survive and GRBs occur 
at z < 1 for HL GRBs and z < 0.05 for LL GRBs, from 
Fig. 9. 

Gamma rays can be produced by not only cosmic-ray 
synchrotron radiation but also neutral pions, charged 
pions, muons and electron-positron pairs generated via 
photomeson production. Gamma rays originating from 
photomeson production have very high energies (e.g., 
> 10 PeV for gamma rays coming from neutral pion de- 
cay) . Such gamma rays cannot avoid attenuation by the 
CIB, CMB and cosmic radio background. In fact, the 
mean free path of 10 PeV photons against pair-creation 
by the CMB photons is ~ 10 kpc, so that we cannot 
expect to detect such gamma rays directly [22:] . Sec- 
ondary electron-positron pairs generated by pair-creation 
are still energetic and up-scatter cosmic background pho- 
tons. These boosted photons can create pairs if they are 
still energetic, and the process repeats itself until the 
energy of degraded photons is in the 1 — 10 TeV range 



[8J, [85J. The mean free path of these regenerated 1 — 10 
TeV photons is longer than 100 Mpc, and they can reach 
the Earth. Whether the detection of these secondary 
gamma rays is possible or not depends strongly on the 
intergalactic magnetic field strength [86j |. As it is large, 
the expected secondary flux becomes highly suppressed 
(see, e.g., |38J and references there in). It is because 
the duration of secondary emission becomes long and the 
emission becomes isotropic (for > 10" 16 G). Only when 
the intergalactic magnetic field is very weak, which can 
be expected in the void region rather than the structured 
region, the secondary emission can be detected. Even 
in such cases, we expect that the secondary flux is sub- 
dominant compared to the primary flux for parameter 
sets in Fig. 9. In fact, the expected fluence originat- 
ing from gamma rays from photomeson production can 
be estimated as £^0 7 ~ 10 ergem -2 for HL GRB and 
-E 2 </> 7 ~ 10~ 8 ergcm~ 2 , respectively (UH^. In addition, 
the secondary fluence from cosmic-ray synchrotron emis- 
sion (estimated by using the same method presented in 
Ref. [38() is smaller than the primary one. Therefore, we 
omit spectra of secondary delayed emission in this pa- 
per [87|. Note that the diffuse gamma-ray background 
from GRBs will be much smaller than the EGRET limit, 
which is expected to be -E 2( & 7 < 10~ 8 GeV cm _2 s _1 sr _1 
[38l [89j|. So far we have considered the internal shock 
model. For other models, see Appendixes D and E. 



V. SUMMARY AND DISCUSSION 

In this paper, we have shown that not only protons 
but also heavy nuclei can be accelerated up to ultra-high 
energies in both of HL GRBs and LL GRBs. We exploit 
the internal shock model, (external) reverse shock model 
and forward shock model. We have also discussed var- 
ious implications for neutrino, gamma-ray and UHECR 
astronomy. Especially, we have studied cosmic-ray accel- 
eration in LL GRBs which may play an important role 
as high-energy cosmic-ray sources. 

Let us summarize this paper below. First, we have 
shown that UHECR nuclei can be produced in both of 
HL GRBs and LL GRBs. 

(Al) In the internal shock model, UHE protons and 
heavier nuclei can be produced in both of HL GRBs 
and LL GRBs. However, the allowed parameter range 
is limited. At smaller radii, most of the UHE nuclei 
are depleted, and a significant fraction of the nonther- 
mal cosmic-ray energy is transferred into neutrinos. On 
the other hand, survival of UHE nuclei is possible only 
at large radii and/or for large Lorentz factors. Typically, 
relatively large radii r > 10 15 cm will be required. For 
HL GRBs with L b ~ io 51 " 52 ergs" 1 , both of protons and 
heavier nuclei can be accelerated up to ultra-high energies 
above 10 20 eV. For LL GRBs with L b ~ 10 46 " 47 ergs" 1 , 
heavy nuclei can be accelerated above 10 20 eV, while pro- 
ton acceleration above such high energies may be diffi- 
cult. In addition, we have shown that the effect of syn- 
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chrotron self-absorption suggested by Wang et al. [TtJ 
is not so important due to the cross section in the non- 
resonant region. 

(A2) In the (external) reverse shock model, UHE nu- 
clei can be produced in both of HL GRBs and LL GRBs. 
Not only heavy nuclei but also protons achieve ultra-high 
energies above ~ 10 20 eV. Whether heavy nuclei can sur- 
vive or not depends on various parameters, and the suffi- 
ciently low photon density is required for survival. That 
is, small E e] , e e , es and n are needed. Conversely, if val- 
ues of these parameters are large enough, in which opti- 
cal/infrared are often expected, survival of UHE nuclei is 
impossible. Note that UHE nuclei produced at internal 
shocks could be photodisintegrated by photons generated 
at the reverse shock, when the ejecta is in the thick ejecta 
regime. Hence, we could not expect survival of UHE nu- 
clei, when we see strong optical/infrared flashes. 

(A3) In the (external) forward shock model, not only 
heavy nuclei but also protons can achieve ultra high en- 
ergies above - 10 20 eV in both of HL GRBs and LL 
GRBs. When the photon density is sufficiently low, sur- 
vival of UHE nuclei is possible. The relevant parameters 
/jv 7 and maximum energy £^ ax depend on time. In the 
ISM case, /at 7 monotonically increases with time before 
the jet break occurs. On the other hand, in the wind- 
medium case, //V7 monotonically decreases with time. 
We have shown that the jet-break effect can be impor- 
tant for survival of UHE nuclei, which was not pointed 
out previously. After the jet break, the photon density 
decreases, so that survival of UHE nuclei becomes easier. 

In the internal shock model and reverse shock model, 
we have to assume that heavy nuclei are contained in rel- 
ativistic outflows of GRBs. It is natural to consider that 
only light elements can be synthesized in such a highly 
relativistic outflows due to the high photon to baryon ra- 
tio [9(J[9l|]. However, there may be possibility that such 
an outflow is contaminated with heavy nuclei that are 
produced by explosive nucleosynthesis [92|, [9^, [94], [95| . 
Also, such heavy nuclei can be entrained from the stel- 
lar surroundings due to Kelvin-Helmholtz instabilities 
and/or oblique shocks 96]. Whether entrained heavy 
nuclei can survive or not is the issue that should be care- 
fully examined. On the other hand, in the forward shock 
model and hypernova model, heavy nuclei will be sup- 
plied by the swept material. They could come from the 
stellar wind of the progenitor star, e.g., the Wolf-Rayet 
star. 

(B) Neutrino astronomy is useful as a direct probe of 
cosmic-ray acceleration, although the detection is not so 
easy. High-energy neutrinos from GRBs have been pre- 
dicted in various contexts. If detected, we can obtain very 
useful information on cosmic-ray acceleration in GRBs. 
In this paper, we have calculated neutrino spectra from 
HL GRBs and LL GRBs in the internal shock model. We 
have shown that the neutrino detection by near-future 
telescopes such as IceCube and KM3Net would be diffi- 
cult, when acceleration and survival of UHE heavy nuclei 
(such as iron) are possible. In the internal shock model, 



for example, we can expect ~ 1 event from one GRB 
at z = 0.1 and ~ 50 events/yr as the cumulative back- 
ground for one of the typical parameter sets (where UHE 
heavy nuclei cannot survive), with large baryon loading 
factors required in the GRB-UHECR hypothesis. On the 
other hand, survival of UHE heavy nuclei leads to much 
smaller neutrino fluxes, as we have demonstrated. Hence, 
second generation neutrino telescopes would be needed 
for neutrino detection in the latter case, although ob- 
servations by IceCube and KM3Net would be useful to 
constrain the GRB-UHECR hypothesis. 

(C) We can also expect high-energy gamma-ray emis- 
sions (including leptonic and hadronic gamma rays) from 
GRBs as well as high-energy neutrinos. We have shown 
that high-energy gamma rays above TeV can escape from 
the source, when UHE heavy nuclei (such as iron) can 
survive. As an example, we have calculated and evalu- 
ated gamma rays coming from accelerated cosmic rays 
in the internal shock model. If the source is nearby, 
we have seen that high-energy gamma rays from cosmic- 
ray synchrotron radiation would be detected by GLAST, 
MAGIC and VERITAS, with large baryon loading fac- 
tors required in the GRB-UHECR hypothesis. (How- 
ever, we should keep in mind that the synchrotron self- 
inverse Compton radiation by accelerated electrons could 
be more important.) Hence, the high-energy gamma-ray 
detection will also be one of the important clues to testing 
the GRB-UHECR hypothesis. When the intrinsic pair- 
creation optical depth becomes high enough (which infers 
UHE nuclei cannot survive) , escape of TeV gamma rays 
becomes impossible and the intrinsic pair-creation cutoff 
should exist. But we may still find signatures of UHECR 
acceleration [78[, although we cannot expect survival of 
UHE heavy nuclei in GRBs. 

In this paper, we have demonstrated spectra of high- 
energy neutrinos and gamma rays, which are expected 
if GRBs are UHECR accelerators, by using detailed nu- 
merical calculations. We have also discussed detectabil- 
ities of such high-energy emission in the near future. 
Not only neutrino and gamma-ray astronomy but also 
UHECR astronomy are becoming important to test the 
GRB-UHECR hypothesis. In Appendix F, we have dis- 
cussed implications of the GRB-UHECR hypothesis for 
UHECR astronomy by exploiting qualitative arguments. 
Future observations of UHECRs may allow us to distin- 
guish between bursting sources and steady sources. If 
UHECR sources are bursting sources such as GRBs, we 
expect to obtain information on the effective EGMF and 
source number density by using the observed anisotropy 
of UHECRs. If we can, we could constrain the local 
bursting rate, for example, which may allow us to dis- 
tinguish between HL GRBs and LL GRBs. 

When we are preparing this work, we found that Wang 
et al. [l7[ studied acceleration of UHE nuclei in GRBs 



independently. In this paper, we have studied in more de- 
tail by making use of the detailed cross section and spec- 
tra of afterglows, and calculate spectra of high-energy 
neutrinos and gamma rays. Furthermore, we have also 



13 



studied cosmic-ray acceleration in LL GRBs. 
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APPENDIX A: COSMIC-RAY ACCELERATION 
AT INTERNAL SHOCKS OF LL GRBS 

In this section, we discuss acceleration and survival of 
UHE nuclei in the internal shock model of LL GRBs. 
We think that it is important to demonstrate whether 
acceleration and survival of UHE nuclei are possible or 



not, because a kind of fine tuning is necessary. In fact, 
the parameter range such that acceleration and survival 
of UHE nuclei is possible (see Fig. 3) is rather limited, 
and we can more easily find parameter sets such that 
acceleration of survival of UHE nuclei is impossible. Note 
that collision radii of LL GRBs arc largely uncertain. We 
can expect r ~ 10 15-16 cm, which implies the time scale 
St ~ 10 2 ~ 3 s HIEE HI. In fact, the observed light curve 
of GRB 060218 is relatively simple and smooth. However, 
the collision radius may be much smaller as suggested in 
Ref. [H. 

As in cases of HL GRBs, we can evaluate the maximum 
energies as (see also Ref. fl4T |) 
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Therefore, we expect that LL GRBs also can produce 
UHECRs, especially that UHE nuclei can be accelerated 
above 10 20 eV. These two conditions (Eqs. (|Aip and 
12[) ) equivalently lead to inequalities (ignoring the term 
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measured in the comoving frame 
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But used parameters are 
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FIG 

L b = lO^ergs- 1 , e" oh = 5 keV, r = 10, r = 9 x 10" cm 
(corresponding to St/(1 + z) — 150 s) and £b(~ es/te) = 1. 
Note that this parameter set is the same one as used in Ref. 
[3] and it implies that a moderate fraction of the energy 
carried by protons goes into neutrinos. 




4 6 

log(s N [GeV]) 

measured in the comoving frame 

FIG 11: The same as Fig. 3. But used parameters are Lt = 
10 47 ergs-\ 4b = 5 keV, r = 5, r = 7 x 10 12 cm, I = 10 11 
cm and £b(~ es/e e ) = 1.1 x 10~ 2 . Note that this parameter 
set is the same one as used in Ref. [64[ and it does not allow 
for accelerating cosmic rays up to ultra-high energies. 
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(1 + z)) as 



1/2 1/2 / r s h(r s h — i) 



1/2 



< 0.55A 4 Z- J 7 rVi 5 rf^ 20 , (A3) 

which is the same as Eq. (2) of Ref. (l4| . As seen 
from the above inequalities, acceleration of protons above 
~ 10 20 eV in LL GRBs will be more difficult than that 
in HL GRBs, because necessary parameters should be 
suitably chosen. For protons to be accelerated above ~ 
10 20 eV, relatively more luminous/magnetized LL-GRBs 
with higher Lorentz factors are required. In fact, GRB 
031203, whose redshift and luminosity are intermediate 
between two dim events (GRB 060218 and GRB 980425) 
and usual HL GRBs, might be such an event [ll|, |97[. 
Eqs. ([Alj) and (|A2|) just mean that UHECR production 
is also possible in relatively lower luminous GRBs such 
as XRFs and LL GRBs, even if LL GRBs do not form a 
distinct population. 

In Fig. 3, we have shown the results in cases where 
acceleration and survival of UHE nuclei are possible. For 
comparison, we show the results for the parameter set 
used in Ref. [14J. In this case, the collision radius is 
smaller, so that we cannot expect survival of UHE nuclei. 
Protons can be accelerated up to ~ 5 x 10 19 eV while Fe 
nuclei up to ~ 10 20 eV (but photodisintegrated). 

In Fig. 5, we show the results for parameters in- 
ferred by Ghisellini et al. [64]. A collision radius is 
small enough, r = 7 x 10 12 cm, so that the photomeson 
production and photodisintcgration efficiencies are rather 
high. Because the photodisintcgration and photomeson 
production processes prohibit acceleration of cosmic rays 
up to very high-energy, we cannot even expect UHECR 
production in this case. 



APPENDIX B: ENERGETICS OF GRBS AND 
UHECRS 

In this section, let us briefly overview the energetics 
of GRBs and UHECRs in order to see that GRBs (and 
hypernovae) could be one of the candidates of UHECRs 
(despite of some caveats). The detailed argu ments are 
found in Rcfs. [3, H H 0, ©, [H B M, WM EqU- In 
the previous section, we have seen that relativistic out- 
flows that make GRBs satisfy various conditions for cos- 
mic rays to be accelerated up to greater than 10 20 eV. 
Here, let us compare the energy generation rate of ob- 
served UHECRs with the gamma-ray energy generation 
rate o f GR Bs. The UHECR generation rate given by 
Ref. [101] from Fly's Eye data is El R dN CR /dE CR w 
0.65 x 10 44 ergs Mpc"V _1 • The UHECR generation rate 
from the recent PAO results [1021 ] is compiled by Der- 
nier HI]. The UHECR generation rate at 10 19 eV is 
E% R dN CR /dE CR w 0.8 xlO 44 ergs Mpc _3 yr -1 - At higher 
energies, the value becomes smaller. 



The radiation energy generation rate of GRBs is given 
by the product of the local GRB rate p(0) (without the 
beaming-correction) and the average (isotropic equiva- 
lent) gamma-ray energy release E™°. First, let us con- 
sider HL GRBs. Recently, Kocevski & Butler [l03| pro- 
vided ^ lkcV ,10MeV] ~ ( 4 T " 7 - 8 ) X 1052 er S S ' in a reSt 

frame bandpass (1 — 10 4 ) keV. Emission at higher en- 
ergies (above the (1 — 10 4 ) keV band) implies that the 
total gamma-ray energy input can be even higher. As- 
suming 13 ~ 2, we can expect E^° ~ (1-2) x 10 53 
ergs. The energy input at high energies is rather un- 
certain, because it depends on the high-energy photon 
index, pair-creation cutoff and possible inverse-Compton 
contribution. Here, we just adopt 2.5 x 10 53 ergs as just a 
reference value for HL GRBs. On the other hand, by us- 
ing logN-logS relationship with the assumption that the 
GRB rate traces, e.g., the star formation rate, various au- 
thors estimated the local GRB rate. In the pre-Swift era, 
the estimated long GRB rate was ~ (0.2 — 1) Gpc~ 3 yr -1 
(see, e.g., [l3[ and references there in). After the launch 
of Swift, many GRBs including high-rcdshift bursts were 
observed, and various authors provided the local GRB 
rate (see, e.g., 0, H[). For example, Ref. || ob- 
tained (0.2 - 0.3)Gpc" 3 yr" 1 assuming that the GRB 
rate traces the star formation rate. However, some au- 
thors claimed that the GRB rate has a faster evolution 
with redshift, which leads to the lower local GRB rate, 
- 0.05Gpc" 3 yr" 1 0,11 II- Next > we sha11 consider LL 
GRBs. However, the radiation energy release and local 
rate of LL GRBs are much more uncertain. The radia- 
tion energy of GRB 060218 is ~ 10 50 ergs, while that of 
GRB 980425 is an order of magnitude smaller. The sug- 
gested local LL GRB rate is ~ 10 2-3 Gpc~ 3 yr _1 which 
are likely to be much higher than the local HL GRB rate 
[Til fl2| . Note that we cannot exclude possibilities that 
there is no such a distinct population, although we as- 
sume the existence of LL GRBs in this paper. Too large 
rates might also be impossible due to constraints from 
obse rvations of SNe Ibc. For example, Soderberg et al. 
104j argued that at most ~ 10 % of SNe Ibc are associ- 



ated with off-beam GRBs based on their late-time radio 
observations of 68 local SNe Ibc. 



Now, we can estimate the UHECR generation rate of 
GRBs from the radiation energy generation rate with 
an unknown baryon loading factor, which is denned by 

tion rate of GRBs can be written as £^ R 



Ucr = £uceU~f. Then, the cosmic-ray energy genera- 
of GRBs 

£acc(47rr 2 ^)r[/ 7 Af. Here, let us introduce R which is de- 
fined as iT 1 = (E* r §^)/S C k, where E* R ^ is the 
cosmic-ray energy input at io 18 - 5-19 - 5 eV. When we as- 
sume p = 2 as the source spectral index of UHECRs, we 
have R = \n(E™ R * / E™ R ) . The minimum energy is not 
well determined theoretically, but we can expect E™ ln ~ 



a fewxTrriArc 



10 1 



eV. We have also assumed that 



cosmic rays can be accelerated up to ultra-high energies, 
so that we can take E™ ax > 10 20 eV. Hence, for HL- 
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GRBs, we obtain 110 
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and for LL GRBs, 
dN C R 
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Therefore, the UHECR energy generation rate is roughly 
comparable to the radiation energy generation rate of 
GRBs (unless, for example, we use the smaller local GRB 
rate inferred from the faster evolution with redshift). 



In other words 

9 dN, 



we can expect 



19 



eV) 



c obdFt( £ ob) for both ofHL GRBs and LL GRBs, and two 
popufations can supply the necessary amount of UHE- 
CRs, when £ acc ~ 10 and p ~ 2. From Eq. (JBTJ, we may 
have the required baryon loading factor £ acc > 20 for HL 
GRBs, which is consistent with Ref. [22|. (If the local 
HL GRB rate is very low, as recently suggested by some 
authors, very large baryon loading factors £a CC ^100 will 
be required, which would be implausible 8]). However, 
its value is not so solid and smaller £ acc could be possi- 
ble, because the local GRB rate (and the total radiation 
energy input and the minimum cosmic-ray energy) is not 
much certain, and we may expect that GRB is one of the 
UHECR candidates at present. 

In the reverse- forward and/or hypernova models, we 
usually use the kinetic energy of ejecta instead of the ra- 
diation energy of the prompt emission. The nonthermal 
cosmic-ray energy is written as £cr. = facc^cj- In the 
case of supernova, e acc > 0.1 is typically needed to ex- 
plain the cosmic-ray flux below the knee. From observa- 
tions of afterglows, we typically have E c j ~ i0 52-53 er g S 
for HL GRBs, while E c] ~ io 50 ~ 51 ergs for LL GRBs. 
In the hypernova model, the total kinetic energy of the 
hypernova ejecta, E C j ~ 5 x 10 52 ergs is expected (see [lH ] 
and references there in). But note that UHECRs can be 
produced only in the high velocity ejecta with T/3 > 0.5, 



which carries E C j ~ z 
models, instead of Eq. 



•^cr. 



dN, 



CR 



x 10 51 ergs [If 
(lB2l). we have 



4.0 x 10 43 ergsMpc J yr^ 1 e a 



dE C R 



For these 



Z?iso 



Pll(0) 



2 x 10 51 ergs / V500Gpc" 3 yr" 
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Hence, UHECRs could be explained in the external shock 
model and/or hypernova model for LL GRBs as well as 
the external shock model for HL GRBs. A significant 



fraction of the kinetic energy must be injected to the non- 
thermal cosmic-ray energy. As in the case of supernovae, 
£acc ~ 0.1 — 1 will be required. Note that the hyper- 
nova model is originally considered in order to explain 
high-energy cosmic rays above the second knee, where 
the cosmic-ray spectrum is formed by the superposition 
of cosmic-rays from hypernova ejecta with various veloc- 
ities. The velocity profile with dE cj /d{Tl3) oc (r/3)~ 3 
leads to the superposed spectral index p ~ 3 (which 
seems too steeper than the typical values of the required 
source spectrum index around the second knee for pro- 
tons, p ~ 2.4 — 2.7). In order to explain UHECRs above 
the ankle only by the hypernova model, we would need 
some additional reason (which was not given in Refs. 
HE E3), such as, e.g., the variable index of the veloc- 
ity profile or the velocity-dependent baryon loading. In- 
stead, the hybrid model (GRBs and hypernovae) might 
be possible. 



The above estimates are based on the assumption that 
the source spectral index is p — 2.0. However, if it 
is steeper, the required energy generation rate becomes 
much larger. In the dip model suggested by @, [Io| . 
the typical spectral index is p ~ 2.4 — 2.7, which is 
steeper than p ~ 2.0 — 2.3 in the ankle model. Such 
a steep spectrum leads to the larger baryon loading 
factor, unless we change Eq^. For example, we have 
i?(10 19 eV) = W(£cr.S ( 1Ql9eV )) ~ 2500 for V = 
2.4 and E™^ ~ 10 115 eV, which can be obtained easily 
from its definition. Furthermore, we have neglected the 
cosmic-ray energy loss at the source. For example, the 
efficient neutrino production leads to depletion of high- 
energy cosmic rays if the shock dissipation at sufficiently 
inner radii. When f pl is not small, the required nonther- 
mal cosmic-ray energy will be raised by ~ 1/(1 — / P7 ). In 
addition, the adiabatic loss and escape of particles should 
be taken into account for more realistic calculations. 



From Eqs. ([Bl]) . (|B2|) and jB3|, we expect that GRBs 
can be the main UHECR sources in both of the internal 
shock model, and the external reverse and forward shock 
model. However, note that too large nonthermal baryon 
loading factors seem implausible because the available en- 
ergy for nonthermal cosmic rays would be limited, e.g., by 
gravitational energy of falling materials of massive stars. 
In additi on, t he high radiative efficiency of the prompt 
emission jl06| might infer that the large baryon loading 
is impossible, as long as the nonthermal baryon energy is 
smaller than the thermal baryon energy. Although these 
possible caveats, current theories cannot answer whether 
efficient UHECR production in GRBs is possible or not. 
Rather, we should test this GRB-hypothesis from obser- 
vations. Current and feature observations of UHECRs, 
neutrinos and gamma rays will give us useful informa- 
tion. For example, observations by AMANDA and Ice- 
Cube have enabled us to constrain the a veraged v alue of 
f PJ under the GRB-UHECR hypothesis fl07U108| . 
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APPENDIX C: THE NEUTRINO BACKGROUND 

It is very important to consider the neutrino back- 
ground from GRBs since it is not easy to detect neutrino 
signals from one GRB event. The cumulative neutrino 
background flux from HL GRBs can be approximately 
evaluated by the following analytical expression using 

r ~ 20 Hlal, 
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(CI) 



c 1 dN' lso 

-min[l,/ P7 ]£ 2 — ^pll(0)/ z 



dE p 

7 x 10- 10 GeVcm- 2 s- 1 str- 1 

f Pll(O) A ( J 
0.05 y V500Gpc- 3 yr _1 



10 



7,50 

(C2) 



where / z is the correction factor for the possible con- 
tribution from high redshift sources, which depends on 
the cosmology. In this paper, we use the standard 
ACDM cosmology with O m = 0.3, fl\ = 0.7; H = 
Tlkms^Mpc" 1 . 

In general, we have to care about other backgrounds 
such as the atmospheric neutrino background when we 
consider the neutrino background. Observations of neu- 
trinos from GRBs have merits that we can take the time- 
and positional-coincidence. Since HL GRBs are bright in 
the gamma-ray energy range, we expect the coincidence 
between neutrinos and gamma rays, which enables us to 
neglect the atmospheric and cosmogenic neutrino back- 
ground essentially. On the other hand, neutrino signals 
from LL-GRBs are very dim in the sense that most of 
the neutrino signals will not correlate with photon sig- 
nals. Only for very nearby bursts, we might be able to 
expect such correlations, and that it requires many-years 
operations. However, we may see neutrino events that 
are positionally correlated with SNe Ic associated with 
LL GRBs. The angular resolution of IceCube for neu- 
trinos is about 1 degree or so, which might be searched 
by the optical-infrared follow-ups with ground-based op- 
tical telescopes [14|. Of course, it is necessary that the 
LL GRB neutrino background makes the dominant con- 
tribution to the neutrino background in this case. 

Note that the GRB-UHECR hypothesis might lead to 
the enhanced cosmogenic neutrino background. Cosmo- 
genic neutrinos are generated when UHECRs propagate 
in the universe, and they are the most promising very 
high-energy neutrino signals. Recent studies have sug- 
gested that GRBs may have the strong evolution with 



redshifts, so that UHECR sources have the strong evolu- 
tion if they are GRBs. In general, the strong evolution 
of UHECR sources leads to the higher cosmogenic neu- 
trino background flux, which may be useful as o ne o f the 
indirect clues to the GRB-UHECR hypothesis (lQ9j . In 
addition, these neutrino signals can give us information 
on the source spectral index, that is, they are useful to 
disti nguis h between the ankle model and the dip model 



APPENDIX D: HIGH-ENERGY NEUTRINOS 
AND GAMMA RAYS IN THE 
REVERSE-FORWARD SHOCK MODEL 

In the reverse-forward shock model as well as the inter- 
nal shock model, neutrino signals are expected. The neu- 
trino emission in the early afterglow phase was studied in 
detail in Refs. dflHO]- We can evaluate the photomeson 
production efficiency as [35[ 



fry ^ 0.088 



L 



ri 6 r- 2 - b 



2 fc ob,10oV 



fa/Elf 1 (E p <E b p ) 
(E p /Elf 1 {El < E p ) 



where e oh is e™ or e^ K or e 



oh 



'ob 



-ob- 



(Dl) 

In the reverse-forward 
shock model, we typically obtain / p7 < 1 except in the 
highest energies. The typical neutrino energy is around ~ 
10 eV rather than ~ 10 15 eV (which is expected for the 
prompt emission). Hence, expected muon event rates by 
IceCube are generally small, and the detection is not so 
easy even by other detectors such as PAO. Nevertheless, 
since EeV neutrinos are produced by ~ 10 20 eV protons, 
the detection of such very high-energy neutrinos is useful 
for diagnosing UHECR acceleration at acceleration sites. 
(For the prompt emission, the typical neutrino energy is 
~PeV, and PeV neutrinos are produced by protons with 
~ 10 17 eV, smaller than ultra high energies.) 

Next, let us consider the forward shock model. From 
Eq. (|D1[) and various quantities derived under the for- 
ward shock model, we can obtain the photomeson pro- 
duction efficiency. By replacing e b and L& with e£ b and 
r 2 (e c £ m )V 2 L £imax , re sp ec tively, we have 



fjyy ~ 2.1 x 10 g_ie ei _ie B / _ 1 S e j,53«o 



(E P /El 



N 0.5 



(D2) 

As seen in the previous section, / P7 increases with 
time in the ISM case (before the jet break), while 
decreases in the wind-medium case. For exam- 
ple, when we assume e acc -E C j = 10 54 ergs and 
Phl(0) = 0.2 Gpc~ 3 yr _1 , the expected diffuse neutrino 
flux at the break energy E b v w 0.05.E* is E 2 V ^ V - 6 x 

10- 11 5 _ 1 e e ,_ 1 e B /2 _ 1 ^e j , 5 3no /2 (/ z /3)GeVcm- 2 s^ 1 sr- 1 . 
In the higher energies, El§ v becomes larger. 

In the cases where UHE nuclei can survive as demon- 
strated in Figs. 6 and 7, we have /Wy ;$ 1 at Ep c ~ 10 20 
eV. From Eq. ([IB]), we obtain f P1 < 4 x 10~ 3 at E p ~ 



17 



10 20 eV. Hence, we can estimate the cumulative back- 
ground neutrino flux as E 2 <& v < 1CP 10 GeVcm _2 s _1 sr _1 
at E v ~ 10 18 eV under the GRB-UHECR hypothesis. 

High-energy gamma-ray emission from the reverse- 
forward shock model is expected as well as the internal 
shock model (see, e.g., [l 12L Ill3] ]). For example, let us 
consider the forward shock model. In deriving Eq. (|19[) . 



let us replace e , a 
(3 = 2 (for p = 2). 
pair-creation at e D b 



= 1 and (3 = 2.2 with e c , a = 1.5 and 
As a result, the optical thickness for 
= 10 12 eVe b.i2 can be evaluated as 



that high-energy gamma rays can escape from the source 
[85| . They demonstrated that not only gamma rays 
from cosmic-ray synchrotron radiation but also sec- 
ondary gamma rays generated via pair-creation by the 
CMB / CIB photons can be detected in the future, if the 
intergalactic magnetic field is weak enough. 



APPENDIX F: GRBS AND UHECR 
ASTRONOMY 



f Nl (E N = 10 20 eV) ' \56 
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where we have assumed e h < e^ h = T 2 (m e c 2 ) /e c oh and 
En > Ejf, which can be realized at t ~ 10 2 s. At later 
time, we can expect e h ~ 10 12 eV > i° h and have 
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f Nj (E N = 10 20 eV) ' \56 
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Hence, we can expect that ~TeV gamma rays escape 
from the source when UHE irons can survive. In fac t, we 
typi cally expect such cases in the afterglow phase [l!2l 
1 1 1 31 ] - Knowing the possible intrinsic pair-creation cutoff 
in the afterglow emission might give us clues to whether 
UHE heavy nuclei can survive or not. However, it may 
be difficult to distinguish the intrinsic pair-creation cutoff 
from other possibilities such as that due to attenuation 
by CMB/CIB photons and that due to the maximum 
accelerated energy of electrons. 



APPENDIX E: HIGH-ENERGY NEUTRINOS 
AND GAMMA RAYS IN THE HYPERNOVA 
MODEL 

In the hypernova model, not only nonthermal photons 
radiated from relativistic electrons but also thermal pho- 
tons (in the optical band) exist at the relatively early 
time. High-energy neutrino emission generated via pho- 
tomeson production is expected, but the detection would 
be difficult Instead, high-energy gamma rays may 

be detected by GLAST and/or MAGIC in the future. 
We can show that the optical depth against these ther- 
mal photons will be smaller than the unity. Although 
photon spectra of hypernovae are very complicated, let 
us make the very simple estimation by using a black-body 
spectrum. We have 



/ 77 (£ob) ~ 0.35 



kT 
leV 
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Therefore, hadronic gamma rays from high-energy neu- 
tral pions, muons and pairs above ~PeV could escape 
from the source. They will be cascaded by CMB/CIB 
photons. After our work was public, Asano & Meszaros 
performed detailed calculations, and they also showed 



The discussions toward UHECR astronomy have been 
recently begun. The direct correlation between UHECR 
sources and observed UHECR arrival directions is one of 
the approaches to study the source properties. This can 
be expected only if the deflection angle due to the galac- 
tic magnetic field (GMF) and EGMF is small enough. 
The latter is poorly known both theoretically and ob- 
servationally. ~ (40 — 95)% of volume within 100 Mpc 
can be regarded as the void region, whose uniform mag- 
netic field is not well-known. The Faraday rotation mea- 
surement of radio signals from distant quasars implies 

S EG dA Mpc ^ 10" 9 GMpc 1/2 [HI EH- In addition, 
there is also the structured EGMF which traces the local 
matter distribution. Observationally, clusters of galaxies 
have strong magnetic fields with (0.1 — 1) ^G at its center. 
This structured EGMF is very important for UHECR 
propagation (see, e.g., [I16l . lll7| and references there in). 
The GMF is relatively well measured, but the magnetic 
field of the galactic halo is not kn own well. This GMF 
also affects trajectories of UHECRs [ll8j . If UHECR tra- 
jectories are deflected only weakly, the deflection angle by 
EGMF can be written as 

V2ZeB EG D 



X 1 / 2 D 



3E N ^D/X 

2.5° ZiEJ^ 2Q Beg. 
ZeB EG D 
2En 

2.6° ZxE N20 Beg,-xiD ioo Mpc- 



100 Mpc 

(Fl) 



(F2) 



Eq. flUD is for the coherent EGMF (9 d D < A) , while 
Eq. ([Fiji is for the EGMF with the coherent length A 
(0 d D > A). Recent PAO results imply the ~ 3° de- 
flection angle, assuming that nearby AGN are UHECR 
sources [Ifl [2l|. The inferred value B EG - 10~ 10 G 
is consistent with the predicted EGMF strength, and it 
could be used to set lower limits on the effective EGMF 
22]. However, in fact, we cannot neglect the effect of the 
GMF. Although the deflection angle due to the GMF 
depends on the model of the galactic magnetic field, it 
can become > 3°. Therefore, lower limits could be set 
only if we see the direction where the effect of the GMF 
is small. In addition, the inferred value is derived from 
the positional correlation between the arrival directions 
of UHECRs and the incomplete AGN catalogue used by 
the PAO collaboration. Therefore, in general, more ob- 
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servations and careful examinations are required in order 
to estimate the effective EGMF from correlation signals. 

Although the GMF can be important for the deflection 
angle, it does not affect the delayed time so much. From 
Eqs. (jFip and (|F2[) . we can estimate the delayed time of 
UHECRs as 



observed at E v — 10 19 ' 9 eV becomes 



Td 



Td 
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EG, -11-^100 Mpc 

Eq. (|F4| is for the coherent EGMF, while Eq. (F3j! is 
for the EGMF with coherent length A. Later we shall 
use the former expression for convenience. Note that 
the EGMF in the above equations expresses the effective 
EGMF field. UHE protons suffer from the stochastic en- 
ergy loss due to photomeson production above the GZK 
energy, so that the dispersion of the arrival time is also 
~ T d - As long as we can use Eq. (IF3[) , we also expect 
multiplets such that the lower energy cosmic-ray event 
precede higher energy cosmic-ra y eve nt, because the time 
delay is statistically distributed jlOOl ]. 

From the delayed time and observed local GRB rate, 
we can estimate the number density of GRB-UHECR 
sources. Let us perform a simple analytic calculation for 
the n umb er density of GRB-UHECR sources. Following 
Ref. [HH, we adopt the top-hat model, where the effect 
of energy losses are approximately negligible when D < 
D c (Ecb.)i and eliminating all cosmic rays coming from 
D > D c (Ecr)- Here, D c (Ecr) is the effective cutoff 
distance. For example, we obtain D c (E p = 10 20 eV) w 50 
Mpc when we define that D c (Ecr) is the radiu s fro m 
which 1/e of cosmic rays can reach us (see, e.g., [l2ll |). 
By assuming that we can observe GRB-UHECR sources 
only during the duration time which is comparable to the 
delayed time ~ Td, we have 



F, 



1 



dN, 



CR 
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AirD 2 T d dE CR ' 
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where Fen is the observed UHECR flux at a given en- 
ergy. Most cosmic rays above the effective cutoff radius 
D c cannot reach us. Hence, the minimum flux -Fq R for 
a given ^ CR is obtained by substituting D c into D. We 
also have the corresponding delayed time rj. The num- 
ber density of GRB-UHECR sources at UHE energies is 
calculated by integrating the differential source number 
density (the number of GRBs per unit observed UHECR 
flux) over the UHECR flux F CR > F£ R . As a result, we 
have n s » (3/5)p(0)r|. Note that, in the GRB-UHECR 
hypothesis, the number density of UHECR sources will 
depend on the observed energy because depends on 
the energy. The number density of GRB-UHECR sources 
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The number density of UHECR sources for arbitrary en- 
ergies can be obtained by calculating D 2 / Eq r . Eqs. (|F6[) 
and (|F7p suggest that determination of the number den- 
sity of UHECR sources and effective EGMF enable us 
to test the GRB-UHECR hypothesis. Acquiring the lat- 
ter is especially not easy, hence one of the ultimate goals 
of the UHECR astronomy. But future observations of 
UHECRs by PAO and TA will give us clues to them via 
observations of correlation signals. 

The source number density can be estimated from the 
study of the small-scale anisotropy. If the observed small- 
scale clustering originates in the same sources, we can es- 
timate th e minimum number density of UHECR sources 
[l22lll23| . For example, recent PAO observations showed 
6 pairs with separation smaller than the correlation an- 
gle scale of 6° among the 27 highest-energy events. This 
leads to a lower limit for the number of sources > 61 
[2~l| . When we adopt D c w 130 Mpc, we can estimate 
the source number density as n s > 3 x 10~ 5 Mpc -3 . 
In more detail, the significance of small-scale cluster- 
ing in the arrival directions of UHECRs can be stud- 
ied by the angular two-point auto-correlation function. 
By simulating the arrival distribution of UHECRs, the 
source number density th at ca n reproduce observational 
results can be estimated [124] . The constraint obtained 
from the small-scale anisotropy observed by AGASA 
is n s > 10~ 5 Mpc~ 3 for uniformly distributed proton 
sources. Despite large uncertainties due to a small num- 
ber of observed events, PAO observati ons in the future 
can lead to the more robust estimation |l25| . 

For example, if we can find that the source number 
density is n s ~ 10~ 4 Mpc -3 , the required effective EGMF 
for a given local GRB rate is £?| G _ 9 Am P c ~ 10 (which 
might be larger the upper limit of the uniform EGMF) 
for HL GRBs and B| G _ 9 Am P c ~ 0.005 for LL GRBs. 
However, since the robust estimation of the source num- 
ber density may be difficult, we may have only the lower 
limit. In such cases, for example, n s > 10~ 5 Mpc -3 
would infer the required EGMF B| G _ 9 Am P c ^ 1 for HL 
GRBs and B| G _ 9 A Mpc > 0.0005 for LL GRBs. There- 
fore, determination of the upper limit on the effective 
EGMF is important in order to test the GRB-UHECR 
hypothesis. If we know that the effective EGMF is too 
weak, possibilities of HL GRBs as observed main UHECR 
sources could be excluded. (However, we could still ex- 
pect to see HL GRBs as rare but very bright UHECR 
sources. In such cases, one HL GRBs could be found 
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as a very bright event among ~ Pll(0)/phl(0) ~ 1000 
LLGRB-UHECR sources.) On the other hand, when the 
effective EGMF is strong enough, HL GRBs rather than 
LL GRBs would be more plausible main UHECR sources. 
It is because the too la rge s ource number density cannot 
explain the anisotropy [125| . Of course, when LL GRBs 
can also produce UHECRs, they could be regarded as 
numerous faint sources, and it is difficult to exclude the 
existence of such faint sources. More generally, we should 
employ the luminosity-weighted distribution rather than 
the luminosity-uniform distribution. Even if LL GRBs 
do not form a distinct population, it would be important 
to know the local GRB rate as a function of E™° for more 
refined investigations. 

Although Eqs. (|F6j) and ljF7|l suggest that future ob- 
servations enable us to test the GRB-UHECR hypoth- 
esis, the above estimation is simple and not quantita- 
tive. In order to present realistic arguments, we need 
the detailed simulation of cosmic-ray propagation, which 
is being planned as our future work. Note that, when 
UHECRs contain heavier nuclei, for heavy nuclei dif- 
fers from that for protons, so that the source number 
density also depends on the composition of UHECRs. 

Before finishing this section, let us briefly describe ob- 
servable features which are expected for GRB-UHECR 
sources. In the GRB-UHECR model, we can expect sev- 
eral prominent features different from those in the steady 
source models. One can be observed in the clustering 
properties. Since the UHECR flux is proportional to D~ A 
as seen in Eq. ()F5[) . the differential source number den- 
sity is proportional to F c ^ 4 as long as we assume the 



luminosity- uniform source distribution [119j . This leads 
to the source number density above i*cR is proportional 

5 I A 

to F CR , which differs from the steady source model 
where we have oc F CR ' . This difference leads to the 
difference in the average number of multiplets, i.e., the 
difference in the ratio of the expe cted number of clusters 
with different multiplicities (123} , In the GRB-UHECR 
model, we may obtain the ratio between the number of 
triplets or higher multiplets and doublets, which is ~ 0.6 
in the GRB-UHECR model, while - 0.33 in the steady 
source model. 

Another prom inent feature is the existence of the crit- 
ical energy |ll9j |. As seen in Eqs. (|F6j) and (|F7|I . the 
source number density decreases with the energy. Since 
D c decreases rapidly above the pion production thresh- 
old energy, will also show the rapid decrease above the 
pion threshold energy. Therefore, there exists the critical 
energy where the observed source number becomes the 
unity. Around this critical energy, we have possibilities 
to see the considerably higher flux than time-averaged 
UHECR flux from all sources. For the source number 
density n s = 10~ 5 Mpc -3 , we expect E c ~ 1.4 x 10 20 eV 
for proton sources. For n s = 10~ 4 Mpc~ 3 , E c ~ 2.5xl0 20 
eV is expected. 

Note that, although we consider the extreme case 
where each population has the typical UHECR luminos- 
ity respectively, there are not necessarily clear differences 
between the two populations. For more detailed stud- 
ies, we should take into account the distribution of the 
UHECR luminosity (i.e., luminosity- weighted model). 
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